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Abstract 

We derive an asymptotic formula a la Liischer for the finite volume 
correction to the pion decay constant: this is expressed as an integral 
over the (3^1^10) amplitude after proper subtraction of the pion pole 
contribution. We analyze the formula numerically at leading and next-to- 
leading order in the chiral expansion. 

1. Introduction The analytical study of finite volume effects is becoming of 
increasing importance as lattice calculations with dynamical fermions approach 
smaller quark masses and aim at higher precision. Since these effects are dom- 
inated by the lightest particles in the spectrum, the pions, and by their long 
distance dynamics, one can study them in the framework of chiral perturbation 
theory (CHPT) [1 . A number of analyses of these effects in different quantities 
have recently appeared in the literature E] • One of these concerned the case 
of the pion mass [2j and has shown that a leading order calculation may receive 
very large corrections from the next-to-leading contribution even for small val- 
ues of the quark masses, whereas even higher order corrections behave according 
to expectations and show a convergent behaviour. This accurate study of the 
convergence of the chiral series has been made possible by the use of Liischer's 
asymptotic formula for the pion mass The formula relates its leading finite- 
volume corrections to an integral over the tttt scattering amplitude in infinite 
volume. Since the latter is known to next-to-next-to-leading order in the chiral 
expansion pj, it is straightforward to evaluate Liischer's formula to the same 
order in the chiral expansion. 



In view of the results for the pion mass, the question arises if one can derive 
similar asymptotic formulae also for other quantities: as we will show in what 
follows, this is the case. In the present article we concentrate on F n , derive an 
asymptotic formula which relates it to the infinite- volume (37r| A^lO) amplitude 
and analyze it numerically using the next-to-leading order calculation of this 
amplitude 0. The results show again large next-to-leading order corrections - 
in this case we cannot explore the chiral expansion further because the two-loop 
calculation of the (37r |v4 M |0) amplitude is not yet available. 

2. The asymptotic formula Denote by F^ t L the pion decay constant in a 
box of size L. The asymptotic formula for AF V = F n ^ — F v can then be written 
as: 

AF * = m^TT / d V e-^^ L N F (zy) + 0(e^ L ) , (1) 
8n 2 M 7T L 

where M> a/ 3 / 2 M w and the amplitude Np{u) is defined as follows. Consider the 
amplitude for creation of three pions out of the vacuum with the axial current: 

(TT 1 ( Pl ) TT 1 ( P2 ) TT 3 (p 3 ) | A J (0) 1 0> = (Pl+P2)»G{S U S 2 ,S 3 ) (2) 

where the superscripts on the pion states and axial current are isospin indices 
and G, H and F are three scalar amplitudes of the variables s±,S2 and S3, with 
si = (pa + P3) 2 an d cyclic permutations 6J. From the amplitude (J2J) one can 
construct the combination which has two of the outcoming pions in an I = 
state (the explicit relation is given below) 

((27r) /=0 7r 3 (p 3 )|^(0)|0) = ( P i+P2),G ( Sl ,s 2 ,s 3 ) (3) 

+ (Pi-P2)»H (s 1 ,s 2 ,s 3 )+p- ilJr F (s 1 ,S2,s 3 ) . 

This amplitude contains a pole in the unphysical region, for (pi + p 2 + p 3 ) 2 = 
Q 2 = M 2 , which needs to be removed before proceeding further. We define 

((27r) /=0 7r 3 (p 3 )l4(0)|0) 5 = ((2vr) /= ovr 3 (p 3 )l4(0)|0) 

iF n T I=0 (s 3 , si — s 2 ) . . 

^ Ml-Q 2 ' {> 

where T I=0 (s,t — u) is the tttt scattering amplitude with isospin zero in the s 
channel. We need the subtracted amplitude in the forward kinematic configura- 
tion, i.e. for pi = —p 2 , S3 = 0, where it becomes a function of one variable only, 
i/ = (s a -Si)/(4M w ): 

^((27r) /=0 vr 3 (p 3 )l4(0)|0) S | pi _ P2 = 2M„uh (v) + M 2 f (u) , (5) 
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where 



h (v) = H (2M v (M n -v),2M ir (M v +u),0) 



and analogously for fo and where the bar on the Fq form factor denotes that it is 
defined after subtraction of the pion pole (the form factor Hq remains unaffected 
by the subtraction). The amplitude Np which enters the asymptotic formula for 
the finite volume corrections to F n is defined as 



The amplitudes H and F can be expressed in terms of F, G and H appearing 
in ©: 

-fo( s l) s 2, S3) = 3-F 12 3 + G23I + ^312 — #231 + #312 , 

#o( s l> s 2, S3) = 3#i23 + - [i*231 ~~ ^312 — G*231 + G312 — #231 — #312] ? (7) 

where Xyk = X(si, Sj, Sk) with X = F,G, H. 

3. Outline of the derivation The derivation of this formula is in large parts 
analogous to the derivation of the formula for the pion mass, due to Luscher [I] . 
In the following we simply outline the necessary steps to prove the formula and 
refer the reader to the paper of Luscher for details. The starting point of the 
analysis is that one can rely on an effective Lagrangian description of the relevant 
physics and analyze these finite volume effects in CHPT. As observed by Luscher, 
the precise form of the effective Lagrangian is never needed in the proof - on the 
other hand, it is very useful to have it available if one wants to understand in 
concrete terms these effects. As was shown by Gasser and Leutwyler one can 
rigorously derive the consequence of chiral symmetry also if the system is closed 
inside a large finite volume with the help of the effective Lagrangian technique 
PQ . In particular the form of the local effective Lagrangian remains unchanged, 
and the only difference with respect to infinite volume calculations comes from 
the propagator for the pion field which becomes periodic in all spatial directions 



where Gq(x) is the propagator in infinite volume. 

The first step in Liischer's proof of the asymptotic formula for the pion mass 
consists in showing that, for a generic loop diagram contributing to the self 
energy of the pion, the dominating finite volume effect is obtained if one takes 
all propagators in infinite volume (G(x) — > Gq(x)) except one, for which only the 



N F (v) = -i (2uh (u) + Mj (u 



)) 



(6) 
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Figure 1: Graphical representation of the asymptotic formula. The wiggly 
(straight) line represents the axial current (pion). The dash on the propagator 
means that it is taken in finite volume (only the contribution with |n| = 1 in the 
sum (JHJ)). Diagram a (b) illustrates the correction to F n (the shift of the pole 
position) due to finite volume. 

terms with |n| = 1 in the sum in (JHJ) should be kept 1 . The sum of all possible 
contributions of this form from all possible loop diagrams gives the leading finite 
volume corrections to the pion mass. The same conclusion is valid also for the 
Feynman diagrams which renormalize the coupling between the axial current and 
the pion - the fact that in this case one of the external legs is the axial current 
instead of a pion does not touch the argument at all. 

The second step in the proof consists in showing, by modifying the integration 
contour in the complex plane, that this leading contribution can be written in a 
very compact form, as an integral over an amplitude (the tttt scattering amplitude 
in the case of the pion mass) defined in Minkowski space, analytically continued 
to complex values of its arguments. Again, the same argument applies also to 
the case of the pion decay constant: in this case, in all possible loop graphs 
that renormalize the pion coupling to the axial current we have to single out 
one internal pion propagator, break it up and put the resulting two external 
legs on shell. The relevant amplitude in this case is the (37r|yl M |0) amplitude, 
as illustrated in Fig. - the weight function which appears in the integral is 
however exactly the same as in the pion mass case. 

The kinematic configuration in which the amplitude must be evaluated is 
also the same and corresponds, for the irn amplitude, to forward scattering. The 
(37r|A M |0) amplitude is however singular for this kinematics because of a pole due 
to one-pion exchange among the axial current and the three outgoing pions. This 
singularity does not belong to the finite volume corrections to and should be 
subtracted. The reason for the presence of this pole can be explained as follows: 
the (7r|yl M |0) amplitude is defined as the residue at the pion pole of a two-point 

More precisely: this concerns only propagators which are contained in at least one loop, 

cf. H 
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function of the axial current and any interpolating field for the pion: 

= J™, - <f)^ ab P{q) (9) 

P{q) = N^q» J dxe^(0\T<Pl(x)Al(0)\0) , 

with N<p the proper normalization factor which depends on the field In finite 
volume both the residue as well as the position of the pole are shifted. Ignoring 
the latter shift corresponds to multiplying Pl{ci) by (Mj — q 2 ) and not by the 
correct {Ml L — q 2 ) and then taking the limit q 2 — *> M 2 . The result, expanded to 
the leading term for asymptotically large volumes, contains a pole for q 2 = M 2 

{Ml - q 2 )P L (q) ~ {Ml - ^Jjff^ = F„ L ~ + . . ■ , (10) 

where AMI = Ml L — Ml is also evaluated to leading order. Since the shift in 
the pion mass is known and given by Luscher's formula, we can subtract the pole 
(which is illustrated in Fig. |IJd) and get the correct finite-volume value of the 
pion decay constant. The result leads to the subtraction prescription given in 
the previous section. 



4. The coupling constant The formula presented here for F n can be 
extended with obvious modifications also to other quantities, e.g. like G n , the 
coupling constant of the pion to the pseudoscalar quark bilinear P % = qi^r l q 



|P*(0)|7r fe ) =5 lk G 7T 



'.ID 



In this case the amplitude that should replace Np{v) in the analogue of Eq. ((TJ) 
is the subtracted P — > 3n amplitude in the limit pi = —p 2 : 



N G {u) 



lim 



<(27r) /= ovr 3 (p 3 )|P 3 (0)|0> 



GJT 



1=0, 



s 3) s l ~ s 2) 



Ml - Q 2 

In this particular case the Ward identity (m = {m u + m d )/2) 

rhG^ , 



F n Ml 



:i2) 



(13) 



which also holds in finite volume, makes the use of such a formula unnecessary: 
from the finite-volume version of Eq. ()13|) one immediately obtains 



AG, 



-'. Rq — Rf + 2R 



M 



(14) 



where Rm is the relative shift for M n . On the other hand, since we have an 
explicit expression for all three relative shifts for large volumes, Eq. (|14|) can be 
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used as a nontrivial check on the asymptotic formulae. Indeed, all three relative 
shifts can be expressed as an integral with the same weight function, and Eq. (JUjl 
can be satisfied only if the same relation holds among the integrands: 

N a (y) = N F (u) - ^-F(u) , (15) 

where F[y) = T /=0 (0, v) is the forward scattering amplitude appearing in Liischer's 
formula for M w . It is easy to verify that this relation follows from the Ward iden- 
tity 2 

-ig^(7r 1 (p 1 )7r 1 (p 2 )7r 3 (p 3 )l4(0)|0) = 7h<7r 1 (p 1 )7r 1 (p a )7r 3 (p 8 )|P s (0)|0> , (16) 

once the limit to the relevant kinematical configuration is taken and if one prop- 
erly accounts for the pole at Q 2 = M 2 present in both amplitudes. 

5. The asymptotic formula in chiral perturbation theory As was shown 
in [2j, the Liischer formula for the pion mass can be used very conveniently in 
combination with the chiral expansion for the irn scattering amplitude. The same 
can be done for F n using the chiral expansion for the infinite- volume (37r|^4 M |0) 
amplitude, which has been calculated up to next-to-leading order in (£]. The 
chiral expansion for the amplitude Np reads 

N F (v) = ^[N 2 F (v)+ZN[(9) + 0(e)} , (17) 

where £ = (M 7r /47rF 7r ) 2 and v = v/M^, and translates into a corresponding 
expansion for AF n 

Rf := ^ = ~ fc/f (A) + eilW + 0(e)] , (18) 

where A = M n L. The integrals I n can be given analytically in terms of a few 
basic integrals: 

= -2B°(X) (19) 
lf(A) = (2l 1 + ^ 2 -3l 4 -05°(A)+(-^ 1 -|l 2 + ^5 2 (A) 

+ ~ (J2°(A) - Rl(X) - WR%(\)) - ^<(A) + |i#(A) + ™R 2 '(\) , 



2 Notice that in the definition of Np, Eqs. (15161) . the (37r|A M |0) amplitude is multiplied with 



^3 and not with as in this Ward identity. 
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where the integrals B and Rf are defined as 



B-(A) = jy«*e-^ = (!) >WA) , (20) 

and 

^ W(A) = {lm fdyy k e-V^^ g (0 {2{1 + m forP^ , (21) 



wi 



ith 3 



g(x) = cr log— — + 2 , tf(x) = - 
a + 1 a; 



2 , a - 1 
— log 1 

ax cr + 1 



(22) 



with a = a/1 — 4/x. These integrals (with the only exception of the primed Rq) 
have already been introduced in [2|. 

We have evaluated numerically these corrections using the following values 
for the chiral low energy constants [Jj: 

4 = -0.4 ±0.6, £ 2 =4.3 ±0.1, £ 4 = 4.4 ±0.2 . (23) 

The results are displayed in Fig.|2]where we plot the modulus of Rp as a function 
of M n for volume sizes between 2 and 4 fm. We have studied the uncertainties 
in Rp which arise from the low energy constants (|2*3*|) and found that they are 
barely visible on the plot - we therefore omit them (in size they are similar to the 
thickness of the lines) . In the figure we compare the evaluation of the asymptotic 
formula to leading and next-to-leading order also to the full one-loop calculation 
of Gasser and Leutwyler jH], which can be given in a very compact form: 

F V>L = F K [1- ^i(A) + 0(£ 2 )] (24) 

where 

, poo 

^(A) = V / dx (™?+™!+"§) A2 , (25) 

Jo 

where the prime indicates that the sum runs over all integer values of Hi, exclud- 
ing the term with all rij = 0. 

In comparison to the pion mass, the finite volume corrections in Eq. ()24j) are 
a factor 4 larger but negative - the sign difference is in accordance with the ob- 
servation that in finite volume chiral symmetry is restored, i.e. the pion becomes 
heavier and its decay constant tends to vanish. Apart from this quantitative 
difference, the numerical analysis gives results which are qualitatively similar to 
those obtained for the pion mass j2J: 

3 The function g(x) is related to the standard J one-loop function through g{x) — 
167r 2 J(xM2). 



6 




M v (GeV) 



Figure 2: The absolute value of the relative finite volume correction Rp = 
Fk,l/F-k — 1 as a function of M v for different volume sizes. We plot the leading 
(LO) and next-to-leading order (NLO) in the chiral expansion of the asymptotic 
formula (j!8|) and also the full one- loop result in CHPT (|24|). The solid lines show 
the sum of the full one-loop result and the NLO correction in the asymptotic 
formula. 



1. the finite volume corrections are exponentially suppressed for large values 
of M n L and become negligible rather quickly; 

2. the leading term in the chiral expansion of the asymptotic formula receives 
large corrections even for the physical values of the quark masses - the 
similarity to the pion mass results makes us however think that the series 
will start to show a convergent behaviour at NNLO; 

3. the leading term in the asymptotic expansion also receives large corrections 
from the subleading ones whenever the finite volume effects are nonnegli- 
gible; 

4. since subleading terms are important both in the chiral as well as in the 
asymptotic expansion, the best estimate of the size of these finite-volume 
corrections is obtained by summing the subleading effects in both expan- 
sions, as shown by the solid curves in Fig. |2j 
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For example, in a recent calculation of F n on the lattice ^1] with dynamical 
fermions a volume of L = 2.5 fm size has been used, and pion masses as low 
as 0.24 GeV. For these values the finite volume corrections evaluated with the 
asymptotic formula to NLO (LO) are 1.5% (1.1%), whereas the full one-loop 
calculation gives 1.6%. Adding both types of subleading effects we find a total 
correction of 2%. In Ref. [T2], L = 1.5 fm and M n = 0.4 GeV were used: in this 
case the full one-loop calculation gives a 3.4% effect, whereas adding the NLO 
chiral corrections we get to 4.5%. For the parameters used in finite- volume 
effects are negligible. 

6. Conclusions We have derived an asymptotic formula for the pion decay 
constant in a finite large volume along the same lines as Liischer's formula for the 
pion mass jlj. The advantage offered by such a formula is a relatively easy access 
to a study of higher order chiral corrections in finite volume effects. We have 
evaluated these numerically and have shown that in F v these corrections are large, 
analogously to what has been found for M v j2] . In the present case we could use 
existing calculations of the relevant infinite-volume amplitude to evaluate next- 
to-leading chiral corrections. Going one order higher in this expansion would 
require the calculation of the (3^1^10) amplitude to two loops in CHPT. 

The asymptotic formula derived here immediately applies (after the necessary 
but obvious modifications) to other similar quantities, like G^. As we have 
explicitly verified, the asymptotic formulae for F n and G n satisfy a Ward identity 
that relates their ratio to M^/m: if one extracts the finite- volume expression for 
M n from this Ward identity one recovers Liischer's formula. The formula applies 
also to the decay constants of heavier mesons, like Fr. In the latter case the 
study of these finite volume effects [9J is of direct phenomenological interest in 
view of the recent application of the lattice calculation of the Fk/F^ ratio to 
the extraction of V us [HI] _ it is worth mentioning that for this application the 
required precision of the lattice result is at the percent level. The same formula 
can also be applied to the decay constants of yet heavier mesons, like fjj or 
fs- In this case, however, the advantage provided by the asymptotic formula 
with respect to a plain one- loop calculation (as recently performed in [Hj) will 
be of practical relevance only if the knowledge of the low energy constants of 
the chiral Lagrangian describing the coupling of heavy mesons to pions ^H] is 
extended beyond leading order. 
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